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ABSIBACT 

Previous  researches  on  the  helix  have  made  simplifying  approxi- 
mations about  the  nature  of  the  "boundary  conditions.   In  this  paper  the 
field  eqxiations  and  houndary  conditions  are  formulated  exactly.  They  are 
then  solved  hy  an  expansion  in  powers  of  the  ratio  of  the  thickness  of  the 
wire  to  the  distance  between  turns. 

The  method  used  consists  of  introducing  a  new  coordinate  system, 
which  is  such  that  a  helical  wire  of  circular  cross-section  is  a  surface  in 
which  one  coordinate  is  constant.  Maxwell's  equations  and  the  electromagnetic 
boundpry  conditions  are  expressed  in  terms  of  this  system.  Since  non-orthogonal 
coordinates  are  involved,  the  equations  cemnot  be  solved  exactly,  but  a  perturba- 
tion procedure  may  be  applied  as  indicated  above. 

The  result  of  the  analysis  is  to  show  that  there  is  a  principal 
mode,  which  propagates  with  the  free  space  velocity  of  li^t  in  the  direction 
of  the  wire.  The  characteristics  of  this  mode  are  studied,  and  they  are  com- 
pared successfully  with  experiment. 


2, 

1,  Introdactlon. 

The  problem  of  the  electronagnetic  modes  of  a  helix  has  been 
treated  "by  making  idealizations  of  the  physical  picture.  Thus,  the  actual 
helical  wire  has  been  idealized  to  a  cylindrical  sheath  conducting  in  a 
fixed  direction  along  its  surface,  and  alternatively  to  a  helical  wire  of 
zero  thickness.  These  idealizations  entail  either  physical  or  mathematical 
difficulties.  In  this  paper  the  problem  of  a  helical  wire  of  finite  cir- 
cular cross-section  is  treated  directly  in  terms  of  a  helical  coordinate 
system.  The  differential  equations  so  obtained  (see  below  eqtiations  3,21) 
are  solved  approximately,  yielding  a  mode  which  propagates  with  the  velocity 
of  light  along  the  direction  of  the  wire.  It  may  be  noted  that  the  coordin- 
ate system  employed  is  not  orthogonal.  The  fields  of  a  helical  sheath  have 
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been  discussed  by,    among  others,  Kompfner     ,  Pierce      ,  Chu  and  Jackson-^, 

Brillouin  ,   and  Phillips  and  Malln-^.      The  fields   of  an  infinitely  thin  helical 
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wire  have  been  discussed  by  Roubine    ,  Parzen  ,   and  Phillips    ,   and  in  an  early 

9 
paper  by  Pochlington  . 

It  has  been  generally  realized  that   the  theory  of  the  helical 
sheath  is  inconsistent  with  the  physical  picture  of  an  actiial  helix.     Thus, 

the  fields  which  are  obtained  depend  on  the  angle  O  ,    in  cylindrical  coordin- 

in0 
ates,   as  e  which  means   that  the  svirfaces  of  consteuit  phase  are  stationary 

in  space  as  we  move  in  the  axial  direction.     However,  the  fields   should  be  ex- 
pected to  turn  at  the   same  rate  as   the  helix.     Further,   the  boundary  conditions 
at   the  sheath  are  satisfied  for  all  Q  ,  while   th^   should  only  be   satisfied  at 
the  wire   itself.      Thus,    the  sheath  corresponds  not  to  a  single  helix  but   to  a 
set  of  adjacent  helices,    filling  up  the  cylindrical   surface. 

The  theory  of  the   infinitely  thin  wire  avoids  this   difficulty  but 
leads  to  another.     To  push  a  finite   current  thru  an  infinitely  thin  wire  requires 
an  Infinite  e,m,f.     Hence  the  tangential  component  of  electric  field,   calculated 
by  means   of  the  retarded  potential,   becomes  infinite  at   the  wire,   while  it  is 
actually  required  to  be  zero  by  the  boundary  conditions.      Since  the  boundary 
conditions  at  the  wire  are  necessary   to  determine   the  propagation  constant  of 
the  fields,    the  theory  is   inherently  inconsistent.     The  difficulty  can  be  avoided 


"by  impressing  different  Taoxmdary  conditions,  Imt  such  a  procedure  is  open  to 
grave  physical  objections. 

The  method  presented  in  this  paper  avoids  these  shortcomings  hy 
obtaining  the  fields  of  a  wire  of  small  hut  finite  thickness.   The  hoiondary 
conditions  correspond  to  perfect  conductivity  at  the  surface  of  the  wire  and 
the  theory  is  entirely  consistent  with  the  physical  picture. 

The  details  of  the  method  are  as  follov/s.  A  special  coordinate 
system  is  introduced,  such  that  the  surface  of  a  helical  v;ire  of  circular  cross- 
section  is  a  surface  on  which  one  of  the  coordinates  is  constant.  Maxwell's 
equations  are  expressed  in  these  coordinates,  and  the  appropriate  ho\indary  con- 
ditions impressed.  The  principal  mode  is  then  calculated  in  detail,  hy  expan- 
sion of  the  fields  in  powers  of  the  ratio  of  the  wire  thickness  to  the  distance 
between  turns.  The  propagation  constant  may  be  found  from  the  boundary  condi- 
tions at  the  wire,  and  the  surface  ciirrent  may  be  determined. 

Once  the  surface  current  has  been  found,  the  fields  away  from  the 
wire  may  be  obtained  by  use  of  retarded  potentials.   The  resulting  formulas 
are  quite  complicated,  but  may  be  separated  into  terms  v-^hich  are  interpreted 
as  circulating  waves,  that  is,  th«  surfaces  of  constant  phase  turn  as  the  waves 

travel.  The  method  of  integration  used  here  is  somewhat  different  from  that  of 

6  8 

Roubine  or  Phillips  . 

We  obtain  the  result  that  thei'e  is  a  principal  mode  which  propagates 
with  the  velocity  of  light  along  the  direction  of  the  wire  under  the  restrictions 
that  the  thickness  of  the  wire  is  small  compared  to  both  the  wavelength  and  the 
separation  betv/een  turns.  The  field  structure  is  given  and  in  particular  the 
VEilue  of  the  longitudinal  field  on  the  axis  is  determined. 

The  method  and  results  are  substantially  those  presented  by  the 
author  at  the  Symposium  on  Traveling  Wave  Tubes  held  in  New  York  on  October  26, 

IS^S   under  the  auspices  of  the  Panel  on  ia.ectron  Tubes  of  the  Research  and  Develop- 
ment Board, 
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2.  Geometry  of  the  Helix 

We  wish  to  determine  a  set  of  coordinates  such  that  a  helical  wire 
of  circular  cross-section  will  he  a  surface  on  which  one  coordinate  is  constant. 
For  this  purpose  we  shall  use  the  following  three  coordinates,  as  shown  in  the 
accompauiying  diagram. 

First,  the  arc  length  along  the  line  of  centers  of  the  wire  cross- 
section  denoted  hy  cK  and  measured  in  such  \inits  that  o<  increased  hy  2n  in  travers- 
ing a  single  turn.   Second,  consider  the  set  of  tuhe  surfaces  concentric  with  the 
wire.  Denote  "by  /*  the  rsidlus  of  the  circular  cross-section  of  a  surface,  measured 
from  the  line  of  centers.  A  tuhe  surface  is  therefore  determined  hy  an  equation 

^=   constant,  and  for  the  surface  of  the  helical  wire  itself  this  constant  is 
equal  to  b,  the  radius  of  the  helical  wire.  The  radius  of  the  cylinder  on  which 
the  line  of  centers  lies  will  he  called  a. 

As  the  third  coordinate  we  choose  an  angle  Q   in  the  cross-sectional 
plane  of  the  tube  surfaces  perpendicular  to  the  line  of  centers  of  the  wire 
measured  from  a  reference  line  determined  as  follows:  The  radius  vector  drawn  from 
the  central  axis  of  the  cylinder  to  a  point  on  the  line  of  centers  is  also  perpen- 
dicular to  the  line  of  centers  and  hence  lies  in  the  cross-sectional  plane.  This 
radius  vector  will  intersect  the  tuhe  surface  in  two  points  A  and  B  which  may  he 
referred  to  as  the  inner  euad  outer  points  respectively.  Now  we  shall  measure 
the  angle  A  from  the  line  connecting  the  center  point  of  the  circular  cross- 
section  of  the  wire  to  the  outer  intersection  point. 

The  three  coordinates  so  defined  completely  specify  a  point  on  any 
tuhe  surface.  These  coordinates  are  single  valued  only  if  the  radius  ^  is 
sufficiently  small.  If  ^   is  increased  to  too  great  a  value,  the  tuhe  surface 
around  one  turn  of  the  line  of  centers  will  intersect  that  around  the  next  turn, 
and  the  angles  ©<  ,  £   are  no  longer  unique.  This  restriction  is  implied  in  all 
subsequent  calculations. 

We  shall  now  derive  the  equations  connecting  the  ordinary  rectangular 
coordinates  x,y,z  with  the  helical  coordinates  ©< ,  /9  i  /*  •  Denote  a  point  on  the 
line  of  centers  by  x  ,  y  ,  z  .  The  pitch  angle  of  the  helix  will  be  called  yj/'    , 

measured  so  that  'lir  =  —  for  a  straight  wire. 


6. 

We  now  hpve: 

X    =  a  cos  oi 

0 

(2.1)  ^o  ^  ^  sinoc 

z     -  a  <\  tan  \/^ 

0 

We  shall  define  three  aaxiliary  vectors,    1,1       and  1      ;    1     will 
"be  the  radius  vector  from  the  central   axis  to  the  line  of  centersT"  By"~itB 
definition,  of  measures   the  amgular  position  of  a  point  on  the  line  of  centers, 
except   that  o<  ranges  from  -  oo   to    oo   .     Hence,   as   far  as  trigonometric  functions 
are  concerned,  oi  =  ^  ,   where  0  is   the  polar  angle   in  cylindrical   coordinates  with 
the  z  axis  along  the  central  axis.      We  thus  have 

(2.2)  i     =  i     cosoc  +  i     8ine< 

We  take  1.  to  he  a  vector  tangential  to  the  line  of  centers.  From 

the  definition  of  a  helix  as  a  line  making  a  constant  pitch  angle  with  the  z- 
axis,  we  have: 

(2.3)  \   =  i<?  cosT^  +  1^  8ln>;^ 

=  -i  8in«C  cos^  +  1  cosoc  cos"«l^  ■♦•  i  sin"*^ 
X  _Z  z 

We  define  a  third  vector  i  hy  the  relation  i  =  i^  x  i  .  We  may 

u  _3i   _*   _£ 

regard  ij.,  ^   ,   an^-  1*  as  a  set  of  moving  axes  with  origin  at  a  point  on  the  line 

of  centers.  Thej  "behave  as  base  vectors  of  a  rectangular  coordinate  system.   The 

vector  i_  is  given  "by  its  definition  as:  , 

(2.U)  i„  =  -i^  since  sin^/^  +  i  cos  «  sin^^  -  i  cos^^ 

The  vector  distance  from  tiie  line  of  centers  to  a  point  on  the  tuhe 
surface  is  given  hy   /3  i^   ,  where  i^  is  a  unit  vector  in  the  direction  from  the 
center  point  to  the  point  on  the  tuhe  surface.     Resolving  i^  into   its  components 
gives,  using  the  definition  of  ^  , 

i^    =  1^"^/  "*■  ^u  ^^^/^ 

B  i     (co8<^    cos^  -  sinot  sin-^   sinV^) 
(2.5)  -^  '^  ^ 

+  1     (8in«   cos^  +  cos  iX  sin^   sinV^) 
-  i     sin/^  cos  V^ 


7. 


The  vector  Joining  the  origin  of  the  tasic  rectangular  coordinate 
system  x,y,z  to  the  point  on  the  wire  is  given  hy 


(2.6) 


(2.7) 


R=i       x+i     y  +  i     z 

-  _x  y  z 

It  also  is   equal  to: 

R=i     X     +i,y     +i     z     +V 

—  xo  /o  zo         r 


since  we  may  resolve  R  into  the  vector  from  the  origin  to  the  line  of  centers 
plus  the  vector  from  the  line  of  centers  to  the  tube  surface.  Equating  these 
two  expressions  and  resolving  into  components  gives: 

X  =  a  cos«<    +  p  (cos  <x   cos^    -  sin  «<  sin^^  sin"^  ) 

(2,g)  y  =  a  sin*t  "♦■  ^(sin  ot  cos^  +  cosoc   sin^  sinV'  ) 

z  =  a  ot  tan\^-  o  sin^  cos"v/^ 

These  are  the  desired  equations  of  transformation. 

3,  Maxwell's  Bquations  in  Helical  Coordinates 

To  transform  Maxwell's  equations  into  the  helical  system  as  described 
above  we  use  the  notation  and  formalism  of  Stratton  .  The  first  step  is  to  cal- 
culate the  unitary  vectors.  They  are  defined  by: 

(3.1)  a,  =  -^—c   r 

^•-'*  '  i    a  u'  — 

Here  u   ,u  ,u-^  denote  *'<  t  (3  %  P    ^'^  *^t   order;  R  =   (x,y, z).     The 
first   derivatives  are  given  by  the  table: 


(3.2)         a)     Xei  =  -  a  8in<3<  -  /O  (sincx   cos^tf  +  cos  «<  sin^  sin  V^) 

y^  =       a  cosoC  +  fi  (cos  oC  cos^  -  sin  oi  sin  ^  sin  V^) 

z  .  =       a  tan^ 

X  =     -  /^  (cos  o<  sin/?    +  sin  oCcos^  sinV^) 

jr^  =     -  /^(sino^  sin^  -  cos  «»^  cos  ^  sini/^) 

z^  =     -  /'cos^  cosV^ 

x^  =     cos  c<  cos  ^  -  sin  <=<  sin/S  sin  V^ 

y^  =     sinc<  cos^  +  cosoC  sinj^  sin^^ 

Zo  -     -  sin^   cos"v^ 
•  J,  A.    Stratton,    "Electromagnetic  Theory",  pp.    38-U7. 


a 
b 
c 
d 
e 
f 

g 
h 


g. 

The  vector  a^  is  given  as  i^  ^  •*■  i     ^  +  i  ^  ,  and  similar 
formula  for  a  and  a  . 

We  now  have  for  the  •unitary  vectors: 

^^•3)     ,  r  ,  -. 

a)     a.      =  -i     |_a  sinot  ^  P  (sin  o(  co^ S  ■*•  cos  o(  sin/?   sin-0-  )] 

+  i   [a  coscK  +  /o(cos  <X  cos^^   -  sin  of  sin^  sin>^  )J 

+  i     a  tan^I^ 
z 

Tj)     a      =  -1    /^(cosot  sin^  +  sin  o<  cos /^  sinV^) 

-i_/«(Bine(  sin/*  -  cosci^  cos^  slnA^) 

-1     />  co%  j3  coba/^ 

c)     a      =  i     (cos  •(  cos/S   -  sin^  sin^  ain'U^) 

"♦■  iy  (sinOt  cos^  +  cosA'  sin^   sinV^) 
-  1       sin/?  cos'\/^ 

It  is   somewhat  more  convenient  to  express  the  transformation  in 
terms   of  the  unit  vectors  of   the  associated  cylindrical  coordinate  system 
rather  than  the  unit  vectors  of  the  rectangular  system.     We  thus  have: 

(3.1^) 

a)     1      =1     cos  ^  -  la     sin  ^ 
X  r  ** 


y  r  ^ 


\>)     i„    =  1,  sin  ^  ■♦•  ifl    cos  O 
(3.5) 


a, 


i)     r*^  =  a*^  +  2a /<3  cos^  *  /*      (cos^  +  sln^  sin'^'V^) 
t)       d^oL  ^  tan-^    ^8in,^sin-t/^ 

When  we  insert   (3.^)  and  (3.5)   into   (3.3)  ^^^  angles  comhine  so  that  the  only  ones 


9. 

appearing  are  d  -  cL  and  ^  .     We  may  simplify  "by  the   fonnulas: 
a)     cos  \0  "  °<)  =     ^— ^=— 


(3.6) 


D)     sin((^-^)=    /-BiM  «i^^ 


The  result  is  that  o(  no  longer  appears  explicitly.  After  these 
simplifications,  v/e  have: 

(3.7)  a)  ai  =  ie    r  +  i^  a  ten ->/^ 

h)  ag  =  -i  -^  sin  ^  (a  ■»•  ^cos^  cos  -/-  ) 

•H^  -^  sin  T^(a  cos  /f  +  />  ) 
-i  ^  cos  jS   COB  -4^ 

c)  a  =  i   i.  /a  cos^  +  ^(cos^^  +  sln^^  sin^V^  )] 

+1^  —  sin  ^  sinA/^-  i  sin  W  cos >/^ 

We  may  find  the  line  elements  without  much  difficulty.     They  are 
given  hy  g. ,  =  a,   •   a     ,     We  have: 

(3.8)  a)     §11  =  r^  -^  a^  tan^\/- 


b 


)     &L2  '^  /^      3in  >/^ 


c)  g,3=  0 

d)  §21  =  /^^  sin^/r 

e)  §22  "  '^ 

f )  §23  =  ° 

g)  §31=  0 

h)  §32  ""  ° 

1)  §33=  1 


10. 


It  will  "be   noted  that  since  g_j^  =  g^g  ~  ^S  *  *?  *^®  ^°^   zero  the  unitary  vectors 

are  not  orthogonal. 

The  unit  vectors  are  therefore  given  by; 

(3.9)  a)  i^  = —  [  i<,  r  +  i^  a  tanV^J 

-    ^rT;^    -      — 

h)  ig  =  -  ij.  ^ip^  (a  +  ^cos^  co8^-^) 

+  1^  Jtl-  (a  cos^  */* )  -  i  cos  A  coB-yf/" 

c)  ^3  =  ij.  7  [a  coB^  +/'(coB^/5  +  sin^  sin^x/^  )] 

♦  i^  r-  a  sin^  sinV'-  i  sini^  cosV^ 
The  volume  factor  g  is  given  hy: 

(3.10)  f=  Ir^  +  a^  tan^V'  ^^  sinl/^  0  j   =  /»  ^  f  r^  +  a^  tan^V'-zO^  sin^^/r} 

=  ^"^  (a  sec'y'+  ^cos^  cos  V^) 


p 

^   sin>^ 


/> 


0 
0  0      1 

We  may  expand  an  arbitrary  vector  F  "by 


(3.11) 


1=  4^  -^^^2*13^ 


g 


The  next   step  is   to  find  the  curl  of  a  vector.     This   is  given  by: 

*     a    riff    -     J^fi) 

Here  f  ,f  ,f  denote  the  covariant  components  of  F.  We  wish  to 

express  them  in  terms  of  the  unit  components  F.  ,  J   ,  T     .   This  is  easily 

•♦  ■'■   -   3 

done,  since  we  have, 

(3.13)                         jj  '^i  i' 
^  '  ^  %j  ^' 

■*  Ibid.,  p.  47,  equation  b3  ' 

*♦  Ibid.  p.  kl,   equation  2k;     p.  kO,   equation  19. 


11. 


Therefore 

(3.1U) 

f. 

1 

= 

We  also 

have : 

(3.15) 

^i 

=/■ 

«ii 

•il 

Putting  these  into   (3.12)  gives: 

(3.16)  V'xF"     A-f^l^l     J-  (. 


s-> 


,2  Vg7. 


a 


9  U^ 


'IL 


fs 


ii 


)^j 


*i^=  f  ^£3 


g. 


ij 


g-i 


au3  VSJj    -    5Z  Vi;/  ''j 


^1^^3^(£. 


g 


2J 


^^     -^J 


au^ 


=ij 


T? 


'JJ 


Formulas  (3. 11)  to  (3.I6)  are  valid  in  any  coordinate  system.  We 
shall  now  specialize  to  the  helical  coordinates  ^  ,  £  ,  P  , 
It  is  convenient  to  define  a  set  of  ODerators 


(3.17) 


a)     D. 


VT*! 


11 


'IJ 


i3i 


g 


2J 


^  ^du^  Vi~      au^  V 


jj 


g 


jj 


-] 


t)     D 


g- 


2J 


11 


g 


VT    '3u^   Vi^ 


^ui 


c)     D. 


.^ia 


go. 


g. 


;1 


11 
g~ 


3J       v^'  leu*    y—- 

In  teruiE  of   the  D's  equation  (3.I6)  hecomes 

I 
J 


(3.18)     V.  F  =  r  K  ^ij    *  ia^sj  *  i3  ^3^]    ^, 


12. 


(3.19) 


The  explicit  forms  of  the  D  operators  are  found  to  he: 

^^■'    f  9     ^31     _     e     ^21   I 


a)     Dii  = 


•11 


r 


11 


^)     D,2= 


/O^  sin  -^/r 


g' 


11 


g 


22 


11    /  9        ^32  _c 


] 


'22 


22 


-   "^ 


dr 


-1^ 


0    D,3  = 


VI" 


/^ 


'33 


9^ 


f2^J 


g 


33 


VI"  ^^ 


d)     D^^ 


vi~ 


g. 


11 


"^  ViT 


p 


e)     D22  = 


'g 


22 


[- 


g 


12 


gl 


•^^22 


P 


VT  ^r 


_S-^  8in^/^ 


13. 


f)     D 


23 


Vg22 


rr  ^f 


k 


3L  ^ 


'g 


33 


ill] 


g)     D 


31 


g 


dec 


Latx 


g 


21 


'hi 


g 


11 
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g 


tv: 


viT         ^'^ 


=  ^f. 
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22 
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c< 


a^ 


g 
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None  of  the  D  operators  involve  o<  explicitly,  "but  only  derivatives 
with  respect  to  oC  .  We  can  accordingly  separate  a  factor  e     from  all  com- 
ponents of  the  electric  and  magnetic  fields,  and  replace  -^ —  ty  i  T  .   This 
is  entirely  to  he  expected,  since  the  fields  should  behave  like  waves  traveling 
along  the  wire.  The  Maxwell  equations  sire: 


(3.2O) 


a)  Vx  B  =  ico/JE 


h)  Vx  H  =  -icoc: 


The  time  dependence  e~  '^       is  implied.     Writine;  B  =  E(  /?    xs  )e 
H  =  '3\AtO}e  ,    the  six  equations  become: 

(3.21)        a)     iu>/<H,     =    iSlL^    .-i^ii^^f    s^    .^.e,*_^eJ 
D)     iayiH^    =    ^[^Vi^    E^      *     ^^slnV-B^    -  irE^J 


c)     iC/^H,     =    ^     [(i£|^.^V^^).^     *^(.r-sin^^)E,l 


d)  -ia>eE«    = 


e)  -icoeB^    = 


dp     °ii   ^       ar 


This  is  a  system  of  six  partial  differential  equations  in  six  unknown 
impressed 
functions.  Since  there  is  no/external  field,  we  can  expect  that  a  field  con- 
figuration which  does  not  vanish  identically  can  only  exist  for  certain  values  of 
the  propagation  constant  A  ,  if  we  assume  the  frequency  coto  be  fixed.  To  each 
characteristic  value  of  V  will  correspond  a  mode  of  transmission. 

In  addition  to  the  differential  equations,  it  is  necessary  to  Imnose 
boundary  conditions  at  the  surface  of  the  helical  wire.  We  shall  assume  the  wire 
to  be  of  infinite  conductivity.  This  implies  that  the  tangential  electric  and 
normal  magnetic  fields  vanish  at  the  surface. 
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Now  the  vectors  i     and  1       are  tangent  to  a  surface  of  constant /p    , 
while  i_     is  normal.      Accordingly  the  "boundary  conditions  are: 

(3.22)  a)     E^     =  0\ 

to)     S^      =  0    >         P  =  h     . 
c)     H^       =  0  J 

The  differential  equations  (3.2l)and  boundary  conditions  (3.22) 
constitute  an  exact  formulation  of  the  prohlem. 

k.   Solution  of  the  Near  Field  Equations. 

The  differential  equations  (3. 21), which  give  the  exact  formalation 
of  the  fields  about  the  helix,  are  far  too  complicated  to  hope  for  an  explicit 
solution.  Accordingly,  we  shall  seek  an  approximate  solution,  based  on  the  con- 
sideration that  the  ratio  of  the  wire  radius  to  the  distance  between  successive 
turns  of  the  helix  is  small.  This  implies  that  the  fields  eo^ound  any  one  turn 
are  not  significantly  affected  by  the  fields  around  adjacent  turns.  The  field 
around  a  section  of  the  wire  may  be  considered  as  composed  of  a  cylindrically 
symmetric  field  produced  by  the  currents  in  this  section  of  the  wire  and  a  dis- 
tortion due  to  the  proximity  of  adjacent  turns.   It  is  the  proximity  effect 
which  we  neglect,  and   as  a  first  approximation  the  fields  may  therefore  be  taken 
as  independent  of  the  angle ^  ,  and  derivatives  with  respect  to<ff  may  be  set 
equal  to  zero. 

To  measure  the  ratio  betv/een  the  wire  radius  and  the  distance  between 
turns  we  may  conveniently  employ  b/a  sec"v|x,  and  then  expand  the  fields  in  powers 
of  this  ratio,  dropping  powers  higher  than  the  first.  For  a  consideration  of  the 
fields  in  the  neighborhood  of  the  wire,  P   is  of  the  same  order  of  magnitude  as  b, 
and  higher  povrers  of  /o  /a  sec  "4^  may  similarly  be  neglected. 

With  this  approximation,  the  quantities  g,,  and  g,  given  by  ^.ga)  and 
(3,10),  become: 

(^.1)         a)-^/i^  =  a  sec>4^ 


b)  V  g   =  a/0  sec  A./^ 
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We  shall  separate  the  fields  into  transverse  electric  and  transverse 
magnetic  modes.   Here  "transverse"  means  having  no  component  along  the  wire. 
We  consider  first  the  transverse  magnetic  case,  and  set  ^  =  0. 

Maxwell's  equations  (3.21)  "become,  with  these  simplifications 


d)  -icocBc,    =  -     i-    ~  /o  H^ 


1    _a_ 


e)-icoeE^     =     -_1        [^^8ln^|r  ^     -  1 T  H^ 


aosec 
We  may  integrate  (U.2a),    olitaining: 


f)  -iwcB^     =     ^^,^^  ^     ll%     H^ 


(^•3)  °    "    f^#^='-*^». 


_   /^sinV/'co8\|/' 


z *«. 


This  restilt,  substituted  into  (U. 2c),  yields  H/,  =  0.   We  now  have 
all  components  of  electric  field  given  in  terms  of  H^  by: 


(U.U)  a)     -icoez^'-    l^    ±-P^6 


>^^       ^^>*:^p-     BJn  >j/ COS  1//       a 


Bp 


c)  -i-.v  ^""r^  H, 
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Maltiplying  U,2b)  "hy  -iCt)€and  substituting  from  (U.1+)  gires: 


2  2 

Here  k     =   copft£ 


We  now  have  on  simplifying: 


9/3  i/ 


/o-h;  *^h^] 


It  follows  that  the  second  term  in  the  "bracket  in  (U,5)  is  of 
order  {/"  /&   seo^)  compared  to  the  first.   Consequently,  it  is  to  "be  neglected, 
since  it  Is  of  the  order  of  magnitude  of  already  neglected  quantities.   With 
this,  we  have: 

/  a 

Equation  (4.9)  is  the  differential  equation  for  the  modified 
Bessel  functions  of  order  one,* 

•  G.  N.  Watson,  "Theory  of  Bessel  Functions",  p.  77. 
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The  ■boundary  conditions  (3.23)  mist  also  "be  simplified.  By  using 
(U.U)  we  see  that  a  sufficient  hoandary  condition  is: 


■2-,_sf- 
2 


Let  us  write  u   «=  °°|  ^  -  k  .  Equation  (U,9)  hecomes: 


a 
(U.ll)       H;'  -^  i  Hy«   -  (u^  +  -i^)  H,  «  0 

This  equation  has  the  two  linearly  independent  solutions: 
(1^.12)  a)     H^     =  kl^iufo  ) 

h)     H^     =  BE^(u^  ) 

Let  us  consider  the  properties  of  these  two  functions,  I,(u^) 
is  a  monotone  increasing  function  of  p   ,   "behaving  like  z)  for^  small,  K^  (uyO  ) 
is  monotone  decreasing,  and  "behaves  like  l/^  for  small.  Hence,  since  we  only 
consider  small  values  of  ^o  ,  we  may  omit  the  solution  AI-|(^/«  )i  which,  when  com- 
paored  to  BK,  (u.^  ),  is  of  the  order  of  magnitude  of  already  neglected  terms. 
Also,  physically,  the  field  must  decrease  as  we  move  away  from  the  wire  which 
again  excludes  the  AI,  (u^  )  solution. 

We  now  must  apply  the  "boundary  condition  (U.IO).  We  have  the  fonmila: 

(U,13)  ^  zK^(z)  =  -S5K^(z) 

from  which  it  follows  that: 

(l+.lU)  0  =  B  u  h  K^(ul)) 

If  the  constant  B  is  not  to  vanish,  u  must  he  a  root  of  equation  (U. lU)„ 

Now  the  principle  "branch  of  K^Cu^b)  has  no  zeros  in  the  finite  portion  of  the  complex 
tt  plane**.  Hence  the  only  pos8i"ble  root  of  (U.lU)  is  u  =  0,  Near  u  =  0,  K^(ub) 
is  approximately  equal  to  -  log  1/2  ut,  so  u  =  0  actually  is  a  root, 
Jrom  the  definition  of  u,  we  accordingly  have: 

(U,15)    T^coB^^    -  k^  «  0 
a^ 

i  =  i  k  a  sec  S^ 


0 


♦Watson,  i"bid,  p,  79 #  equation  5. 
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These  roots  are  associated  vd.th  waves  traveling  with  equal  velocities 
in  opposite  directions.   It  will  now  "be  showi  that  they  actually  propagate  with 
the  velocity  of  light  along  the  wire.   This  is  easily  seen.   The  fields  contain 
as  a  factor  exp   -i(&i  t  -  Poi)   which  is  equal  to  exp  -lcu(t-  ^  ^^°  f  ol).  cC   ia- 
creases  hy  2rT  in  going  the  distaince  along  the  wire  2n   a  sec  a^  ,  whence  we  may 
write,  calling  s  the  actual  arc  length  along  the  wire,  exp  -iC4;(t  -  s/c).  This 
is  a  wave  going  at  the  velocity  of  light  along  the  wire. 

The  next  step  is  to  determine  the  actual  field  components  which  may 
"be  most  conveniently  expressed  in  terms  of  the  total  current  carried  "by  the  wire. 
The  surface  current  density  is  given  "by: 

(u.i6)  J  =  -  a  ^  5  • 

Here  n  is   the  tinit  vector  normal  to  the  wire,   for  which  we  may  teifce 
i_   .     Hence: 

(U.I7)  J  =   (ig  X  i^)  H^ 

The  vector  (i-  x  i   )  is  also  of  unit  magnitude  and  is  equal  to: 

/>.  lax           t     ^  *       -     4     /^sin/?   sin'vA'cos'vX^     .    .       cos  \^  ,     ^  \  ^  ^        ^      , 

(1+,18)  ^23  r     — r —  ^    — r^  ^^  + /ocos^  )  +  i^  sin)^ 

To  the  degree  of  approximation  to  which  we  are  working,  this  "becomes: 
(^.19)  ^o^  ^1     "     ^0    cos-»^+  i^   sin>/^ 

The  total  current  is  obtained  "by  integrating  (U. I7)  around  the  wire, 
or  what   is  the  same  thing,   since  H^  is  Independent  of  ^  multiplying  "by  2TT"b.     Hence: 

(U.20)  I     =   2n"b     BK^(uh) 

-*  ~     anh      Kj^Cuh; 

^/«       2rrh        ICj^Cut; 

If  u/*  is  sufficiently  .small,  we  may  replace  K,  (u x?  )  "by  l/u/^  ,  and 
K,  (u"b)  "by  l/u"b.  Then: 

(U.21)  H^  =  ^ 

*   Stratton,  loc.  cit. ,  p.  37» 
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This  expression  is  the  same  as  that  for  the  field  around  a  straight 
wire*,  as  is  to  he  expected,  since  we  have  neglected  the  proximity  effect  due 
to  the  coiling  of  the  wire.  The  other  field  components  are  now  easily  found 
from  {k.k)   and  (4.13)  as: 

(4.22)        ,  ^  II       ^o^^/^  ^  1  I  2,  ,        1 

*^  ^^  =  -  UJT  2Rt     K^(uh)        ^  iZJi    2Fb     ^^  1°S  2     ^/* 

^^     ^=1^    2Fh         g^(uh)       r-  Bi^V-go»^—  i|[7^2^  ^^^g  |u^-f  sm^cos^A- 

°''    */o"    ~  cue    2nh    Kj^(uh)   '^  '¥6  2n/» 

From  these  equations  we  see  that  1!^  and  £^  are   negligihle  off  the 
wire,  since  n  is  zero.  Actually,  "because  the  conductivity  is  finite,  u  is  not 
precisely  zero,  and  hence  we  have  retained  the  forms  (4.22)  a,  "b. 

This  completes  the  determination  of  the  propagation  constant  of  the 
principal  mode  hy  use  of  the  near  fields.  We  shall  next  use  the  current  density 
(U.I7)  to  find  the  far  fields  hy  use  of  retarded  potentials. 

5.  The  Fields  off  the  Wire 

To  find  the  fields  off  the  wire,  we  shall  use  the  retarded  potentleds, 
in  particular,  the  retarded  Hertz  Tector  given  "by:** 


(5.1)        ^  i     f  e^'^-^" 

Eence  we  have  written  (fi  -  R«|  for[(x  -  x')^  *(/  -  y')^+(z  -  k')^J 
If  we  suhstitute  from  (2.8)  we  have: 


1/2 


(5.2)  /e  -  R«f  =  r^  +  a^  -  2a^  cos  {Q  -^)  +  (z-a  tan^fc')^  */^ 

-  2p*1(t  cos(e  -«')-a|cos^  *+^T  sin(0-«n' )8ini^-(z-a  taax/ic' )coa^ij)^'j 

Here  H'  =   i  oc* ,  £*  t  ^' )  represents  a  point  on  the  surface  of  the  helix. 
To  the  degree  of  approximation  to  vdiich  we  are  working,  we  may  drop  terms  Involving 

•  Stratton,   Ihid,  p.    535,  Equation  52. 
••Stratton,    ibid,  p.   U3I,  Equation  U4. 
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/»,  and.   obtain: 
(5.3)  (fi  -  R'l  ^  =  r^  +  a^  -  2arco9  (<?-«')  +  (z  -  a  tan1//p<')^ 
The  current  vector  J  is  given  'by: 

(5.1+)    J  =  (i^cooV*V  '^^^^  2Sb  ^^^'^ 

In  order  to  perform  the  integrations,  we  wish  to  eliminate  the 
dependence  of  the  vectors  i  ^'  smd  i  ,  on  the  coordinate  o<  .  This  may  be  done 
as  follows: 

(5.3)  1..  «  i. 

1  ,  »  -1  sin  ^  +  1  cos  ff 

9  ^       y 

=  -(i     cos  5  -  i^sind  )  sin  or'  +  (i  8in^+  i^  cos  ^  )fc08  ex' 

=     i     sin(^-*<0*i^    cos   (<?-P<') 

The  volume  element  becomes  the  surface  element  2rrb  a  sec  W  d.<^  ,  whence 
the  integral  reduces  to: 

(5.6)  [[=     YmuT    ^  *  secl/^  Ado^Ti^  sin  {6 -ci)  co8V^+  i^coa  i0-J)cos>l/ •*■  i  6in)^e    %-R'f 

exp  ik  I r^  +  a^  -  2a  r  cos   (^-«^)  +  (z  -  a  tan  1// (X  ')^J 

We  shall  make  the  change  of  variables    of  -    — — -j-  =of,  write 

actan  I//  • 

/    =  ka  sec  1^  ,  and  write^-    — r — zrr    -  4^  •     Tiie  integral  becomes: 

^5' 7)  ^  lkasec|f|&(+|^] 

11=  ^_L_.    I  a  sec\J//  dPC  ij.sin((/>  -  o<)co6V/   +i^cos(0-«)  008^^+  ^g^inyU       /jK-R'l 


-  ~  imue 


•o» 


1/2 


exp     ik   jr^+  a^  -  2a  r  cos   (^  -^)  +  a^  tan^yc<^J 

I  ir  z 

Let  us  %rriteT  =  k  csc^/",  whence  the  integral  becomes   e  times  a 

function  of  <J>  above.     How  the  surfaces  on  which  <^  and  r  are  constant  are  helices. 
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whence  we  see  that  we  have  a  representation  in  terms  of  circulating  waves, 
turning  at  the  same  rate  as  the  helix.  • 

The  integral  majr  he  evaluated  "by  use  of  the  Integral  representa- 
tion of  Somaerfeld**  ^j 

We  insert  this  into  (,'5*7 )»   fi^d  obtain: 

(5.9)  7r=       i — I  a  sec-v^  e*"^  */doC  (i  sin(<* -o<lcos-^  +  i  ^  cos(</>-o<)cosl^+  i  sin-v^)e' 
—      HTTCcpe  -«•     •-  r         "^  "  8 

/, -/TrrT  /.       A       -^Tv-?    a  tan>/^K/ 

CX  vr  +a  -2a  r  co8(^ -eiiO     e  ~ 


J 


XdX    J 


0  '^ -1^ 

The  X   integral  is  over  a  path  indented  helow  at  \  =  k.   This  is  equiva- 
lent to  assuming  that  k  has  a  small  positive  imaginary  part,  which  is  set  equal  to 
zero  after  the  integration.   The  double  integral  converges  absolutely)   and  we  may 
change  the  order  of  the  integrals.   The  Bessel  function,  J  ,  n«y  be  expanded  by 
means  of  the  addition  theorem*** 

(5.10)  J  (xVr^  +  a^-2a  r  cos(</> -«)  /  =  Z7  J  (>vr)  J  (Xa)e^''^*^  ""^ 
0  -OS  n     n 

The  series  that  results  converges  uniformly,  and  we  may  integrate 
term  by  term.  We  obtain  after  these  transformations: 

^Icxfi^  sin((D-«)cosV*  i^  co8(0-^)cosV'*  ,^3,^^1  in(^ -«HlkasecV'^-V>^dtanV^K| 

The  OC  integral  may  be  evaluated  by  elementary  means.   It  is  convenient 
to  replace  sines  and  cosines  by  exponentials.   Calling  the  integral  U  ,  there 

•  S.  A.  Schelkanoff,  "Electromagnetic  Waves",  D.  van  Nostrand  Co.,  I9U3,  p.  I4O9 
**Stratton,  loc.  cit. ,  p.  576,  Equation  17;  Watson  loc.  cit. ,p,Ul6, Equation  U. 
***Stratton,  loc.  cit.  ,  p.  373,  Equation  8. 


23. 


results: 


(5.12)  ^=iiA\t     iS^Zt/^e^C^^DC^^-*^)  ..l(n-l)(0-<>t)) 
—  -ao     /_r  2x   ^  / 

i^  £082^^^i(n+l)(^  -Of)  ^^l(n-l)(^-<v) 


,ik  a  secV^  -  Vx^-k^  a  tan-»^/«/ 


fora: 

(5.13)     / 


All  the  integrals  may  be  eraluated  as  special  cases  of  the  general 


^iuo^    .T/o</^^^  ^ 


.2^2 


(5.U) 


Ihis  gives  upon  substitution  into  (5.12) 


.v5a? 


U  =Vx'-k' 


.ant[l,22f±('. 


l(n+l)^ 


l(ii-l)i* 


2  ^2^    2^      2 


(X  -k'^)aSanV*(ka  sec^A-n-l)2       (x2_k^2^^^2^^(j,^33,^.^,^^ 


J 


■♦•  ig  cosv/ — 5 — - — 5^ 

—  MX  -k^)a^anV+ 


<^ 


c(n-/}0 


V'-H(ka  seo/'-n-l)'^  (X^-k2)aSanV*(ka  sec^^+D^ 


) 


■♦■12  sin-vl/- 


in0 


(X^-k^)a^tanV'*'(ka  seo/'-n)^ 
2  .     2 


] 


Let  us  factor  out  a     tan  >^  from  the  denominator  and  write 
^n  °  la  tan>^     ~  ^  '^8C>')     -  k       .     We  may  note  that    p^  in  our  notation  is   the  p 
of  other  writers  since  k  C8c>/r  =r.     With  this  shorthand,  we  have  finally: 
2.-2!  .   /       i{n+ll<»  i(tt-l)^ 

i(irH)^ 


(5.15) 


~         a  tan 


*'»     ""'Vl    2     2- 


~~  X'^+    p^    J 


2U. 


Substituting  this  value  for  U  Into  (5.II)  glvew: 
(5.16)  o.      ^ 


[^^(- 


'O 


l(n+l)^  ^i(n-l)^ 


^  Xi-1  ^  *  Pn-1 


; 


+  1        cos-fl-^S __     +     _| ^ ) 


In*     7 

n 

The  various  integrals  appearing  here  may  all  be  evaluated  as  special 
cases  of  the  general  formula*: 

(5.17)    f'xd^J^UrjJ^Ua) 


2    2         °  ^n^^'^  ^n  ^^^'^    '^  "^  * 
X  +  p 

I^(pa)  K^  (pr)    r  >  a 

2  2 

If  p     has  the  negative  vsdue  -q  ,  we  replace  I   (pr)  K  (pa)  by 

l/?^iJ  (q.r)H'   (q.a).     We  shall  only  write  down  results  for  r  <  a,   since  those 

for  r  >  a  can  be  obtained  by  interchanging  r  euid  a  in  the  fornalas  for  r  <  a. 
Thus  we  have: 

•  Watson,   loc,   cit.   p.    U^29,  Equation  5. 
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This  foumila  may  "be  simplified  by  writing  n+1  «=  n'ln  the  first  terms  of 

the  quantities  la  parentheses ;  n-1  ■  n'  In  the  second  terms,  and  then  dropping 
the  primes.   !Illhe  result  Is: 

The  convergence  of  series  of  this  type  has  "been  discussed  "by  Phillips  . 

It  Hiay  easily  "be  shown  that  If  r  f'  a,  the  terms  of  the  series  decrease  geometrically 

if  n  Is  sufficiently  large.  For  r  =  a,  the  series  converges  conditionally  if 

&f   z/a  tanV".  In  the  case  r  =  a,  0=   s/a  tan-y^  ,  the  series  diverges.  However, 

this  corresponds  to  the  neighborhood  of  the  helical  wire  Itself,  and  the  eqpproxima- 

tlons  under  which  (5.19)  waa  derived  are  not  valid  in  this  region  and  the  divergence 

does  not  affect  our  discussion. 

The  next  step  is  to  o"btain  the  field  quantities.  These  are  found  from 

* 
the  Herts  vector  "by  the  operations 

(5.20)  a)  B  »  V  (V  •!)  f  k^I 

•b)  H»  -iccfc  VxTT 

The  fornolas  for  divergence  and  curl  in  cylindrical  coordinates  are: 
'     V    U  r     9r      "  r       r     dJ  qz   "z 

^i.(?|f  'iT*-^^ir,) 


•  Stratton,   loc.   cit.,  p.   U3O 
••Ibid,  p.   51. 


26. 


The  formulas  may  "be  simplified  ty  use  of  the  recurrence  relations 
for  Bessel  fonctions.     The  formula: 

(5.22)  l|-rypr)«     p(li(pr)*    ^^ 

enables  us  to  write  out  V  *  li  as: 

/-nl     ,(pr)  nl       (p  r)  -j 

*  -V-  IP.  (      Tl  Vl  V)  *      °*p/  Vl(Pn»)} 

Combining  the  terms  in  curly  brackets  gives  for  the  entire  bracketed 
e^qpression: 

(n-H)l   (p  r)  -J 

-Oi*lV^  *  i7^-^)^nn  (Pn^>] 

*  2  sinV^  l(r  -  --^^)l^(p^r)  Z^(p^a)] 


The  expression  in  c\irly  brackets  here  may  be  simplified  by  the  re- 
curring  relations: 


(5.25)      a)   I^(«)  *yi^(«)  -  I^.iU) 


♦  Watson,  loc,  cit.,  p.  79,  Equations  3  and  U. 


27. 
When,  these  are  inserted  we  ohtain  for  the  curly  bracketed  expression: 

(5.26)  I^^r)  [K^.i(P^a).K^^^(p^a)} 
Another  recurrence  relation** 

(5.27)  K  ,(p  a)  -  Z  ^-(p  a)  =  -  -^  K  (p  a) 
^■^*  '  '  n-1  -^^n     n+1  *^n       p  a    n  -^n 

■^n 

reduces  (5.26)  to; 

■^n 
We  now  have  as  the  slsipllfied  value  of  (5.2U) 

(5.29,      1^  .  s^  ,^  .  ^^  r^i^^r)  t^(^j.y 

*  2  slnlt-  l(r  -  j-^)  I^(p_^r)  I^Cp^aJ 

The  first  term  ceoicels  the  tern  — r-^ — r  In  the  factor  JT-  — r r  » 

a  taJi>K  a  tanV 

so  the  final  result  is: 

(5.30)  2  slnV-   ir   I^(p^r)  K^(p^a) 

(5.31)V-F=  ^^l  cscV-  e^^*_r   e  a  tan^       .   ^^^  i^(p^OK>^a) 

Instead  of  finding  all  field  components  at  all  points  of  space,  we 
shall  only  compute  X  (r=0),  which  is  the  relevant  quantity  for  heam  interaction. 
By  the  prescription  of  Eqxiations  (5.20),  we  have: 


_9_  V  ,n  ^^2JJ 
This  is  equal  to: 


(5.32)       \-Tl      '^"^^    z 


[ 


i^'^-r^^)  •  2ik*k2.  2  3^^J 


a  tan-^f' 


**  Itid,  p.  79,  Squation  1. 


2S. 

At  r  =  0,  Ijj(Pj^r)  is  zero  except  for  n  =  0,  so  the  series  reduces  to 
the  single  term: 

(5.3U)  V°^  '  Tmbe  i  "^V^e^""'  ^^(Po*)  2  (k^sinV--  kr  ) 

-  -  l^^I  •  2k2  cotV  e^^«  ^•''^  I^(ka  cot  >^  )  . 

6.  ConclttBion 

Ve  shall  now  compare  the  results  of  this  theory  with  those  of  pre-vioue 
theories.  Firrt,  we  have  ohtained  a  node  ^ich  propagates  with  the  phase  velocity 
C  along  the  wire  and  the  phase  velocity  c  sin^Aalong  the  axis.  This  has  also 
"been  ohtained  in  previous  theories.  However,  the  restrictions  on  the  validity 
of  the  calculations  are  different.  Here  the  restriction  is  that  h/a  sec -v^ should 
he  small.   In  the  helical  sheath  theory  this  value  of  phase  velocity  is  only 
obtained  for  ka  1  surge,  *  In  the  zero  thickness  wire  theory  a  completely  artificial 
boundso-y  condition  is  required.  Also,  we  have  shown  that  only  two  modes  exist, 
propagating  with  eqtial  velooitles  in  opposite  directions.  The  sheath  theory  lea. 4 
to  an  infinitude  of  modes,  of  which  only  two  are  propagated.  We  have  not  herein 
considered  attenuated  modes,  which  are  not  geherated  under  our  restrictions. 

!Ehe  longitudinal  field  along  the  axis  is  given  "by  equation  (5.3*^). 
For  the  purpose  of  calculating  the  interaction  of  the  fields  with  an  electron  beam, 
this  is  the  relevant  quantity.  It  may  be  used  as  a  factor  in  Pierce's  parameter  C, 
which  determines  the  gain  of  a  traveling  wave  tube. 

In  our  calculations  we  have  obtained  the  propagation  constant  correct 
to  lowest  order  terms  in  b/a  sec-\^.  Now  it  is  a  consequence  of  the  general  theory 
of  guided  wave  propagation**  that  the  propagation  constant  of  a  mode  is  stationary 
with  respect  to  small  variations  from  the  true  field.  Hence,  the  linear  tenn  in 
the  series  expansion  of  T  in  powers  of  b/a  sec -x^  will  vanish,  and  we  may  expect 

*Cha  and  Jackson,  loc.  cit.,  Equation  3^  ff> 
♦♦Schelkunoff,  loc.  cit.,  p.  3SU,  J.  Schwinger,  unpublished  notes. 


29. 
that  a  more  accarate  calculatioxi  vlll  yield  a  r  of  the  form: 
(6.1)  JT  =  k  c8cV^  [l  +  Cj^("b/a  sec  V' )    J 

The  numerical  coefficient  c^  will  pro'bably  be  on  the  order  of  mag- 

11 
nitude  of  jr  to  ^    . 

The  major  results  we  have  obtained  are  the  propagation  constant 
and  the  presence  of  only  one  mode.     These  have  been  confirmed  by  experiment. 
Because   of  the  difficulty  of  calculating  the  fields  by  use  of  the  circoilating 
wave  expansion.   Section  ^,  no  attempt  has  been  made  to  compare  the  field  strength 
within  the  helix  with  experiment. 

He  may  sum  up  as  follows.      The  method  herein  presented  obtains  all 
the  results  of  previous  theories  of  a  cold  helix.     Moreover,   it  obtains  them  in  a 
consistent  manner  without   employing  any  artificial  boundary  conditions,   but 
directly  from  Maxwell's  equations  and  the  standard  perfect  conductivity  boundary 
conditions. 

As  far  as  mathematical  interest   is   concerned,    the  major     point  is 
the  use  of  a  non-orthogonal  coordinate  system.     Despite  the  no n- separability  of 
the  differential  equations,   it   is   still  possible  to  obtain  all  desired  results. 
The  method  may  be  extended  to  treat  the  higher  approximations  and  obtain  the 
attenuated  modes. 
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